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' We study the spectral properties of a chiral random banded matrix (chRBM) with elements de- 

CIh' caying as a power- law Tiij ^ This model is equivalent to a chiral ID Anderson Hamiltonian 

D , with long range power-law hopping. In the weak disorder limit we obtain explicit nonperturbative 

CO ■ analytical results for the density of states (DoS) and the two-level correlation function (TLCF) by 

mapping the chRBM onto a nonlinear a model. We also put forward, by exploiting the relation 
between the chRBM at a = 1 and a generalized chiral random matrix model, an exact expression for 
the above correlation functions. We give compelling analytical and numerical evidence that for this 
value the chRBM reproduces all the features of an Anderson transition. Finally we discuss possible 
C , applications of our results to quantum chromodynamics (QCD). 

^ • PACS numbers: 72.15.Rn, 71.30.+h, 05.45.Df, 05.40.-a 
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I. INTRODUCTION 



The spectral properties of a disordered system, namely, a noninteracting particle in a random potential, are strongly 
. affected by the global symmetries of the Hamiltonian. In the ergodic/metallic regime, the spectral correlations of 
^ ' generic quantum complex systems with translational invariant spectrum fall in three universality classes ('standard' 
from now on) attending to time reversal and spin symmetries. The simplest representative of each universality 
class corresponds to an ensemble of matrices with the appropriate symmetry and random Gaussian entries. Typical 
properties of the spectral fluctuations of these matrices include level repulsion and spectral rigidity. In the literature 
they are usually referred as Wigner-Dyson (WD) statistics [ij. 

In recent years new ('nonstandard' from now on) universality classes have been put forward by relaxing the transla- 
tional invariance condition y,|2l,yiS|- These classes are related to additional symmetries of the system. For instance, 
1^ ' for (random) matrices with block structure 

. where C is a general nx {n + 1/) real, complex or quaternionic matrix, the eigenvalues come in pairs of ie^ or are zero. 
' This discrete symmetry is usually called chiral and induces an additional level repulsion around zero which results in 
^ , different spectral correlations for eigenvalues near zero (the origin) and away from zero (the bulk). In the bulk the 
I spectral correlations are not affected by the block structure and WD statistics applies. 
' ^ , The motivation to study these nonstandard symmetry classes comes from different branches of theoretical physics. 
In the context of QCD, the Dirac operator in a chiral basis has a similar block structure. It turns out that in the 
infrared (ergodic) limit, the eigenvalue correlations of this operator do not depend on the dynamics details of the QCD 
Lagrangian but only on the global (chiral) symmetries of the QCD partition function j^. Thus random matrices with 
. 5^ ' the correct chiral symmetry of QCD (termed chiral random matrices) |^ accurately describe the spectral properties of 
I the QCD Dirac operator up to some scale known as the Thouless energy 0. As in the standard case, depending on the 
invariance under time reversal symmetry and spin, there exist three chiral random matrix ensembles. In the context 
I of QCD each universality class is related to both the number of colors and the fermionic representation considered. 
Another interesting problem that falls beyond the standard classification is that of the spectral correlations of a 
mesoscopic disordered system composed of both normal conducting and superconducting parts. In this case, in an 
appropriate basis, the mean field Bogouliobov-de Gennes Hamiltonian can be expressed as 

where A is a matrix representing the pairing field and h is a matrix representation of the disordered free Hamiltonian. 

The Hamiltonian Ti. can be effectively modeled as a random matrix provided that the phase shift due to reflection 
in the NS interphase vanishes on average In this situation the gap at the chemical potential disappears and thus 
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pseudo particle excitations appear at arbitrarily low energies. The resulting spectral properties also fall apart from 
WD statistics. In this case four new universality classes are found attending to spin and time reversal symmetries. 
Other systems with similar nonstandard symmetries include random flux models (here disorder is placed not in 
the site but in the link (bond disorder)) and bosons in random media 

As a general rule (valid for all standard and nonstandard symmetry classes), localization effects tend to gradually 
erase the impact of symmetries. However, the transition to localization is not the same in both cases. In systems with 
nonstandard symmetries the localization properties are in general sensitive to details as the microscopic form of the 
random potential ^10] or the number (odd or even) of lattice points considered In certain situations 12, 13, 141, 
for energies close to zero, the DoS diverges and eigenvectors remain delocalized even in the strong disorder regime 
in ID or 2D. By contrast, in the standard case the spectral correlations are insensitive to the microscopic details of 
the dynamics provided that the disordered potential is short-range. The eigenstates of a free particle in a disordered 
medium in less than 3D are localized in the thermodynamic limit for any amount of disorder. In three and higher 
dimensions there exists a metal insulator transition (Anderson transition) for a critical amount of disorder. 

For finite systems, the dimensionless conductance g — Ecj ^ {Ec = h^jtc is the Thouless energy, tc is the classical 
time to cross the sample diffusively, and A is the quantum mean level spacing) is a useful magnitude to quantify 
the deviations from WD statistics due to wavefunction localization. We recall that WD statistics applies in the 
metallic/delocalized regime g ^ oo. Nonperturbative corrections due to a finite g ^ 1 (weak disorder) were recently 
evaluated by Andreev and Altshuler llSl by mapping the localization problem onto a supersymmetric nonlinear a 
model O. They managed to express the TLCF in terms of the spectral determinant of the classical diffusion 
operator. As disorder strength increases, localization becomes dominant g ^ 1 and eventually the system undergoes a 
metal insulator transition. Unfortunately the above analytical tools cease to be applicable in this region. Numerical 
simulations of 3D short range Anderson models suggest that, at the Anderson transition, the wavefunction moments 
Pq present anomalous scaling with respect to the sample size L, Pq = J d'^r\^{r)\'^'^ (x L~^''^'^~^\ where Dq is a set 
of exponents describing the transition |17|. Wavefunctions with such a nontrivial scaling are said to be multifractal 
(for a review see Ref.jl^). Spectral fluctuations at the Anderson transition (commonly referred as 'critical statistics' 
[l9|') are intermediate between WD and Poisson statistics. Typical features include: scale invariant spectrum 
level repulsion, and sub-Poisson number variance i21||. Different generalized random matrix model (gRMM) have been 
successfully employed to describe critical statistics |22l l23l | . 

The study of these gRMM have shown that critical statistics and multifractal wavefunctions can be reproduced in 
the weak coupling regime <? 3> 1 by allowing long range hopping in the original Anderson model. Actually power-law 
long range hopping, far from being a mathematical curiosity, appears in a broad range of systems: glasses with strong 
dipole interactions in real space [2J| , the evolution operator of Floquet systems |2^ with a nonanalytical potential or 
in Hamiltonians leading to classical anomalous diffusion [2^ . Although the introduction of long range hopping dates 
back to the famous Anderson's paper on localization these models did not attract too much attention until the 
numerical work of Oono |2^|23 and the renormalization group (RG) treatment of Levitov in the context of glassy 
systems. The main conclusion of these works was that power-law hopping may induce an Anderson transition in any 
dimension provided that the exponent of the hopping decay matches the dimension of the space. 

A ID version of this problem, a REM with a power-law decay, was discussed in Ref. |23| . By mapping the problem 
onto a nonlinear a model with a nonlocal interaction, it was shown that for a 1/r band decay the eigenstates are 
multifractal and the spectral correlations resemble the ones at the Anderson transition. In this case the TLCF in the 
g ^ \ limit is expressed in terms of the spectral determinant of a classical anomalous diffusion operator related to 
classical ballistic diffusion and 1// noise [s^ ■ Recent investigations [sj have further corroborated the close relation 
between this random banded matrix model and the Anderson model at the metal insulator transition. 

In this work we want to study the interplay between long range disorder and chiral symmetry. Our main aim is to 
provide a detailed analytical account on how the transition to localization occurs in systems with chiral symmetry. 

The organization of the paper is as follows. In Sec|nl we introduce a chRBM with power-law decay (equivalent to 
a ID Anderson model with power-law hopping and chiral symmetry). After mapping it onto a nonlinear a model we 
investigate the localization properties by using the renormalization group formalism. In Sec lllll analytical results for 
both the DoS and the TLCF are obtained in different domains. We show for a special value of the power-law exponent 
{a — 1) the resulting expression is greatly simplified. In Sec lIVI by exploiting the relation between the chRBM for 
a = 1 and an exactly solvable random matrix model, we give an exact result for the DoS and the TLCF valid in the 
weak disorder regime. It is shown that the spectral correlations for this special value resemble those of a disordered 
system at the Anderson transition. We also conjecture an expression for the DoS and the TLCF for arbitrary power- 
law exponent. These analytical findings are supported by numerical results from direct diagonalization of the chRBM. 
Finally, in Sec0 we discuss applications of our model in the context of QCD 0- 
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II. THE CHIRAL RBM WITH POWER-LAW DISORDER 

We study the spectral properties of an ensemble of Hermitian chRBM given by 

), (3) 

where wi,2 are nx n Hermitian matrices. This form is related to Eq.(^ at = by a unitary transformation. In this 
paper we restrict ourselves to the v = case. The matrix elements (aji,2)y are independently distributed complex 
Gaussian variables with zero mean and variance 

{\{uj,^2h\')=a\r)= ] (4) 

1 + (r/ro) 

where rp and a are real parameters and r = \i — j\ is the spectral distance. The choice of complex matrix elements 
corresponds to a matrix model with unitary symmetry. Due to the chiral symmetry, the eigenvalues of Eq.Q come 
in pairs of ±6^. This feature induces thus an additional level repulsion around zero and consequently the spectral 
correlations for eigenvalues near zero (the origin) and away from zero (the bulk) are essentially different. At the bulk 
the spectral correlations are not affected by the block structure and coincide with the nonchiral version of Eq.lPJ 
which has been intensively studied in recent years |23l l3l| . 

In this section we study the localization properties at the origin by mapping the chRBM onto a supersymmetric 
nonlinear a model The above chRBM can also be interpreted as a ID chiral Anderson model with long range 
disorder. In this context the parameter tq measures the amount of off-diagonal disorder. Thus the tq 3> 1 {ro <^ 1) 
limit corresponds with the weak (strong) disorder regime. In this paper, due to technical reasons (large-rg allows us to 
use the saddle-point approximation in the derivation of the nonlinear a model) , we are focused on the weak disorder 
limit. Concerning the power-law exponent a, we are mainly interested in the range 1/2 < a < 3/2. It is known that 
the spectral correlations for the nonchiral RBM case a < 1/2 are well described by WD statistics We expect this 
holds also in the present chiral case and we do not consider it in this paper. On the other hand, as shown below, for 
a > 3/2 our model is similar to a standard ID short range chiral Anderson model recently investigated in Ref . |3^ . 

A. Supersymmetric nonlinear a model 

We use the supersymmetry method to derive the sigma model. We shall follow Efetov's notations and conven- 
tions. As usual, the first step is to express the correlation function as a ratio of two determinant. For the one point 
retarded Green function. 



^ ' e+-H 2 dJ ^ ^ 



(5) 

,7=0 



where = e + iS is the energy with small imaginary part and the generating function Z[J] is given by 

_ det(e+-7^ + J) 

The source field Jij is a 2n x 2n complex matrix. The above ratio can be written as a Gaussian integral over a 
supervector with both fermionic and bosonic components as 

Z[J] = [ V{^/j, V^) exp [ii:{e+ --H + fc J)V'] , (7) 



where k = diag(l, —1) in superspace, -0 is a 4n-component supervector and "0 — "0^. 

Following the standard method, we perform the spectral averaging over Ti. and introduce the 4x4 supermatrix Q 
using the Hubbard-Stratonovich transformation. After these steps are carried out, the generating function is given by 



Z[J]) = J PQexp iA-^)^]StTQ^QJ + str In (e+S^ + kJY., + iy^Q 



(8) 



where Aij = a^(|i — j\), Aq = n ^X]"j=i^y' Pauli matrix in chiral space and the Q-matrix verifies 

{Qj'Sx} = 0. In the limit Aq ^ tq 3> 1 the integral over Q can be performed by the saddle-point method. It is 
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not hard to show that in this case the solution of the saddle point equation is Qi = at e = 0. Consequently the 
saddle point manifold = 1 is parametrized as Q = TYi^T where T is the matrix with symmetry specified below 
and T — T^^. In addition to that, the massive mode changing the saddle point should also be integrated out. Unlike 
the standard case, the effect of these modes plays an important role in our case. The matrix Q is thus parametrized 
as 



= T,(E, + 5Q,)%, SQ, = Sq,^,, (9) 

where 6qi is a 2x2 supermatrix. In order to proceed we first perform an expansion in powers of 6Q up to second order, 
and then integrate over it (it is just a Gaussian integration). After these manipulations we find (Z[J]) — J exp(— i^) 
with 

F=^J2 ^^J^*"" Q'Qj " i ^ ^''""^^ f,T,j:,stT f,Tj^, + £1 ^ str S,Q, + ^ str J^jki:,Q,, (10) 

ij ij i ij 

1-1^ 



where Rij — {A^^)ij — n^^Sij J2kii^^^)kh and A — 7rV^Ao/2n is the inverse of the DoS (mean level spacing) at e = 0. 

We remark that the above results can be applied to any a(r) provided that Aq is large enough. For technical details 
we refer to Ref . |32j| where a a(r) with exponential decay was considered. Compared with the sigma model for the 
nonchiral RBM, we found an additional contribution having the double-supertrace form (second term of Ea. (|l(J|l '). 
This term was derived by integrating over the massive modes. We recall that in its derivation it was crucial that the 
massive modes were parametrized as in Eq.ljHJ. 

In order to proceed we have to make a gradient expansion of the kinetic term in Ea. (|10|l . As in the standard case 
[23|. this expansion depends on the value of a. For a > 3/2 

^^i?.,strQ,Q, ~ _^^(j_j-)2^,^.str(VQ.)' 



(2) 

dxstr (VQ(x))^ (11) 



where 



oo 



4Ao 



1 



l + ir/ro)^" sin: 

2 



2a 



1 + (r/ro)2" 

4(2) „• JR_ 

-l-=rl-^- (12) 
Ao sm ^ 

We see that the integral in AJ, is well-defined for a > 3/2. After a similar expansion in the double trace term above, 
we find the normal diffusive a model for a > 3/2 

j dxsir T,^Q{x) + dxsti (VQix))^ ~ ^ J dx (str Q{x)VQ{x)j:^f , (13) 



2AV 



where for the sake of simplicity we have neglected the source term, and V — n is the 'system volume'. The coupling 
constants b and c are given by 



b Ao ""^sin^' c AAl 8 ^sin^' 

ZOL ^' 2a. loL 



(14) 



For 1/2 < a < 3/2 the gradient expansion requires a special care since Aq is not well-defined. Taking into account 
higher order terms in the gradient expansion (as in B,ef. |23j| V we find the following expression for 

.2 , T ^ 1-COSgr _ „ 2a-l,„,2a-l 



r ^'^TXrX^^^/ dr - ~arg°-^|gr-\ (15) 

' l + ir/roY" Jo l + (r/roy°' 
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where Cq = 2 Jp°° dx{l 
given by 

F 



where 



cosa;)/x^" is a numerical constant. With this substitution Ea. H10() for 1/2 < a < 3/2 is 

- ^ I dxstrE,Q(x) + l I g|gp"-istrg(q)g(-g) 

~h J ^^^l9l'"~'strQ(p + q)g(-p)S,strg(p'-9)0(-y)S., (16) 



^0 



T 1 sin 1 

2a ^ 



^ „2q 
-"a ' 

"2ir 



2a 



(17) 



In this case, unHke for a > 3/2, the nonlinear a mode (I16|l corresponds with a process of anomalous diffusion. However, 
we note that the relation 1/b = iAo/c ~ ro/c holds in both cases and 1/c <C 1/b. Thus the double trace term (last 
term of F) will not affect the spectral properties to leading order in l/rg. Its role in the localization properties of the 
chRBM is discussed below. 



B. Renormalization group equations 

By using the above nonlinear a model Eas. p3|l and p6|l . one can discuss the localization properties of our chRBM 
in the thermodynamic limit. This can be done by investigating the running of the coupling constants b and c under 
the RG flow. We compute the RG equations to leading order in b and b^/c. For the normal diffusive case, Ea. (|13|l . 
the parameters b and 6^/c involve the following propagators 



U{q, uj) = — , n2(g, LJ) = q^n\q, Lj) (For a > 3/2), (18) 



1 

q'^ — luj 

respectively. On the other hand, in the anomalous diffusive case, Ea. (|16|l . the propagators are given by 

n('7:'^) = ^^r4 -' Tl2{q,Lo)^\q\^''-'nHq,co) (For 1/2 < a < 3/2). (19) 

By observing the momentum integrations of the propagators 



— --^n(g,tj) - — In/, / — — n2(g,w 
A/; (Stt)" 27r J^/; (Stt)" 



^^In/, (20) 



where A is a cutoff and Z '--^ n > 1, we conclude that for normal diffusion the logarithmic dimension is d = 2. For 
anomalous diffusion, since our chRBM is in essence one dimensional, the logarithmic dimension corresponds to the 
case a — 1 (generally, d = 2a — 1). We note that it is a straightforward task to extend the derived cr model for the 
chRBM to arbitrary dimension. 

In the normal diffusive case, the perturbative calculation of the RG equations was first done in Ref. [12 by the 
method of the e-expansion. The RG equations are given by 

^ db , ^ dc 1,^16^ ,„ , 

f3b = -— = eb, p^ = - =ec+—c\ ( ^ , (21) 

d/i d/i OTT ZTT C 

where fi = hil is the renormalization scale, /3f,_c are the beta functions for b and c and C is the wavefunction renormal- 
ization. The above equations are immediately applied to the anomalous case by replacing e — d— 2 by e = 1 — {2a — 1). 
The RG equations are easily solved for arbitrary e We see that, in the domain e > 0, the system is within the 
extended phase. For e < 0, b and 6^/c goes to infinity and consequently a transition to localization is expected. 
However we note that in this limit the perturbative expansion breaks down since b, c increase under the RG flow. 
For e = 0, 6 remains constant, c goes slowly (~ 1/log/) to zero, and 6^/c goes to infinity. That means that at the 
origin and for exponentially large volumes the eigenstates should be delocalized and the DoS is expected to diverge. 
Of course this divergence is somewhere cut off before the origin since random matrix theory (which should be valid 
in the deep infrared limit) predicts a vanishing of the DoS at the origin. 

Finally we remark that the above pseudo-divergence is beyond the current numerical capabilities since b <^ c and 
1 /c ~ log I ^ log n for n ^ 1 . 
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III. LEVEL CORRELATION FUNCTIONS 



In this section we compute the DoS and the TLCF close to the origin in the 5 S> 1 hmit by using the nonhnear a 
model introduced in the previous section. As mentioned previously, in systems with chiral symmetry there is a clear 
distinction between the spectral correlations near zero (the origin) and away from zero (the bulk). In the bulk the 
spectral correlations are not affected by the block structure and coincide with the nonchiral version of Eq. Q which 
has been intensively studied in recent years |23l l3ll | . After reviewing the result in the bulk, we examine the origin for 
the anomalous diffusion case 1/2 < a < 3/2. For normal diffusion, we refer to Ref.|22j. The case a = 1, related to 
the Anderson transition, is worked out in detail. 



A. Review of results at the bulk 



The use of the supersymmetry method permits an analytical evaluation of both spectral properties and eigenfunction 
statistics |23 in a certain region of parameters. In the thermodynamic limit it was found that the eigenfunctions 
are multifractal for a = 1 and localized (delocalized) for a > 1 (a < 1) |2^. The spectral correlations in the limit 
g — ttEc/A ^ 1 {Ec the Thouless energy) can be expressed through the spectral determinant of a classical diffusion 
operator |l5i |. For the unitary ensemble 

1 1 92 , , , cos(2z)„9, , , , 

" 52(2^„)2(2a-l)+^2' (23) 

n=l ^ ^ ' 

where R[z) — ^ l^p(t\) pit-i)) — 1 is the TLCF, p(e) is the DoS at energy e, A is the mean level spacing, 1/2 < a < 3/2, 
and the energy z is expressed in units of A as z = 7r(ei — e2)/A. Although this result was derived assuming 0^1, 
in Ref.js^ it was shown that it is valid for all z. 

We mention that the spectral determinant above corresponds with a process of anomalous diffusion '33| with 
(|r|) oc xhe asymptotic behavior of i?(2) ~ ^-2+i/(2a-i) _^ jg power-law and the dimensionless 

conductance g ~ n^"^" increases (decreases) with the system size n for a < 1 (a > 1). Both the scaling of 5 and the 
TLCF resemble those of a weakly disordered conductor in d = 2/(2a — 1) dimensions '23|. For the special case a = 1, 
V[z) = (z/2g)sinh"\z/25) and R{z) can be explicitly evaluated as 

^W = 1-^('-#!7tV- (24) 



z^ \sinh(z/2(7) 

This correlation function reproduces typical features of the spectral correlations at the Anderson transition as level 
repulsion for z <C 1 but sub-Poisson number variance ^^(i) ~ xL for L ^ 1 (see Sec lIVI for a definition of the number 
variance). 

Although the RBM is defined in ID, generalizations to d dimensions are straightforward. In that case the product 
over n runs over all possible n — {ni, . . . ,nd). As mentioned previously, according to Levitov's results [2J|, the 
properties of these power-law hopping models are similar in different dimensions provided that the dimension of the 
space is equal to the power-law decay. It is tempting to guess that a power-law decay a in d = 1 is similar to a decay 
da in d dimensions. 



B. Results at the origin 



After reviewing the properties at the bulk we move to the spectral correlations close to the origin. In the limit 



g = ttE^/A 



^2-2q 



/b 



„2a- 



-2a 



> 1, we calculate the DoS and the TLCF of the chRBM Q by using the 



nonlinear a model as given by Ea. (|16|l without the double trace term. As mentioned previously, since b/c ^ l/?'o ^ 1; 
the double-trace term can be safely neglected. 
The DoS is given by the expression 



(Pie)) 



1 



4:AV 



Re / VQ 



dxstr kYizQ{x) 



(25) 
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where F is defined in Ea. H16() . The perturbative calculation corresponds to expand the Q-matrix as 



Q{x) = 



1 + iP 
1-iP 



P = 



t 
t 



t = 



p ib 



(26) 



where a,b are real variables and a,p Grassmann variables. The resulting expansion is in powers of the (anomalous) 
diffusion propagator 



n(g,e) 



1 



27ri:)|g|2"-i 

where D = An/nb. Finite g corrections to the DoS are given by 



(P(^)> 



1 

A 



1 



(27) 



(28) 



Since the expansion is not well-defined for the zero mode {q — 0), this contribution should be removed from the 
above expression and treated in a separate way [s^l- We remark, as was noticed in Ref. J2J, the mean level spacing 
A corresponding to the ergodic regime is modified by a finite g as 



1 

A 



1 



Re / VQ 



1 



dxstr kYizQix) 



'F[Q] 



1 

A 



iRe(^n(q,, 



g#0 



(29) 



where Q denotes the nonzero modes. This quantity is relevant for us since we are interested in level correlations of 
unfolded variables. Thus we define the DoS scaled in terms of the renormalized mean level spacing 



pi(z) = A(p(e = Az/7r)). 



(30) 



We note that, as shown in Sec llVI the renormalization of A is important to find agreement between the results of this 
section and those of Ref. [s^ . 

The zero mode contribution is treated nonperturbatively and gives the ergodic result 5 — s- oo |36| 



pi{z) p[ \z) 



'-iJ^{z) + jUz)). 



(31) 



As in the nonchiral case a proper parametrization of the Q-matrix is required to derive this expression. It can 
be incorporated into our Q-matrix as 



Q{x) = TQ{x)f, 

where the supermatrix T (T is the inverse of T) parametrizes the zero mode as |37 

T = UTatJ, 



(32) 



Tn = 



U 



cos ^ 
— i sin ■! 



u 
u 



1 2 



u — exp 



Of 
i9B 



c 

TJ 



(33) 



where — tt < Op < n, < 9b < 00, and ^ and 77 are Grassmann variables. Using this parametrization, we integrate the 
zero mode first and, then, treat the nonzero modes perturbatively. This was done in Ref. (3^ for the case of normal 
diffusion. In the present case, the zero mode part is unchanged since the kinetic term, second (and third) term in the 
free energy H16() . does not include the zero mode. For the nonzero mode, the anomalous diffusion propagator Ea. H27() 
is used for the perturbative expansion. Thus after integration of the zero mode, we are left with. 



1 7T d 

(p(e)) - x + TTX^^^l— / dt{t^zo){Jo{t + zoAB)Jo(t-zoAF)-Jiit + zoAB)Ji{t-zoAF)) 
A 2A dzo J^^ 



kin ' 



/kin 



1 

2V 



dxstr ■ 



' Zo 

l±k 



S,[Q(a;)-S,], 



= / (• • •) exp 



1 



dq 



-71 7T^\<ir-'stvQ{q)Qi-q) 



46 / 2tt 



(34) 
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where zq = ne/A and 1/A is the purely perturbative contribution 1)29(1 . We note that this expression is vahd up to 
second order in 1/g. The ergodic result can be easily recovered by putting Ap s = 0. 

We are now ready to evaluate Eg . ((34() in different domains. In all cases the limit 5 3> f is assumed. For z g, 
known as the Kravtsov-Mirlin (KM) domain 34] , the Bessel functions can be expanded in powers of zA ^ z/g as 



J{t - zA) 



I- zA 



dt 



2 df^ 



Jit), 



(35) 



to obtain 



-Re 



Changing the variable from zo = 7re/A to z = 7re/A, we finally are left with 



Piiz) 



1 



0.a I ^ d 2 d 



P?\zh 



where p^^\z) is defined in Ea. H31|) and Cq, is the momentum summation for periodic boundary conditions, 

1 



(27rn 



l4a-2 ■ 



(36) 



(37) 



(38) 



We now move to the Andreev-Altshuler (A A) domain J^] z ^ 1. In this limit we cannot expand the Bessel 
functions. Instead we use their asymptotic form, 



Mz) 



1 



dxe^ 



Noting that Ji(z) = —Jq{z) we find 



Piiz) 
Fk{z) 



Re 



1 



1 

8^ 



8z^ 2z 







1 




j cos z + 




• sin z 






-8l + - 





dq 



46 / 27r' 



|2a-l 



SivQ(q)Q{-q) 



iz 
2V 



dxsti k'EziQ — S^). 



(39) 

(40) 
(41) 



We remark that the presence of the supermatrix k = diag(l, —1) breaks the supersymmetry in Fk{z). Keeping terms 
up to second order in the P-matrix we obtain 



cos 2z , 



1 



(42) 



where the spectral determinant I'(z) is given by Ea. (|23|l . We note that the ergodic limit corresponds to put V = 1 
which corresponds with the asymptotic form of the exact result Eg. pifl . 

We now consider the TLCF. In this case the analytical calculation is along the lines of the DoS though more 
technically involved. Here we present a summary of results. For technical details we refer to Ref.|32| where the case 
of normal diffusion was discussed. We start by defining 

P2{zi,Z2) = A2((p(ei = Azi/7r)p(e2 = Azs/tt))), (43) 

where ((p(ei)p(e2))) = (p(ei)p(e2)) — {p{^i)) {p{^2)) is the connected part of the TLCF. In the KM's domain zi_2 <C g, 



^2(^1,2:2) 



452 



d 



d 



dzi dz2 



1 + -r^ Zi- V Z2-^ — + — ^ z^—^ + 2ziZ2-^-^ h Z2-^ 



8.9 



92 



d d 



52 



9z2 



■ dzi dz2 



'■dzl 



K{zi,Z2) 



where 



K{zi,Z2) = ^/^^ ^1 {ziJi{zi)Jq{z2) - Z2Jo{zi)Ji{z2))- 



(44) 



(45) 
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The ergodic limit is given by ^2(^1, ^2) = —K'^izi, z-^). In the AA's domain z\^2 ^ 1, 

P2(zi,z2) - iRe^(nUnl) + ^Piim^(n+ + n_) + ^2?2im^(n+-n_) 

[X»ir'2(2?+2?r^-l)cos2(zi+Z2)+X'iI?2(2'^2?+^-l)cos2(2i-Z2)] 



8ziZ2 



^ [1 + I?i2?22?+2?I^ cos 2(zi + Z2)] - 5^ [1 - 2?i2?22?-2?+^ cos 2(zi - Z2)] 



2(zi + Z2)2 ' --^-+-- — ■ -^J 2(Z1-Z2)2 

+ ^-^-^ (X>i sin 2zi - X»2 sin 2z2 ) , (46) 

^1 ~ ^2 

where I?i,2 = 2?(^i.2), = V{{zr ± Z2)/2), Eq.(|lSIl, and ^^,2 = n(g,ei,2), n± = n(q, (ei ± e2)/2), Eq.(E7I). We 
remark that the first term of the r.h.s. is obtained by a purely perturbativc calculation. 

We stress that in the the common domain 1 ^ z ^ g the results in both regions AA and KM coincide. It is 
worthwhile to note that in the unitary limit z, zi + Z2 00 our results also coincide with those of the nonchiral 
version of our model: pi{z) — > 1 and 

i?(^r,^.) - 1 + - 1 + 5: - -JL^ + £|^(fl^P^ (47) 

Pl{zi)pi(Z2) 2 2{zi-Z2Y 2{zi-Z2Y 

Noting n(g, e/2;g) = 2n(q, e; 2^) and I?(z/2;g) = 'D{z;2g), we see that this is the AA's result (|22) for the unitary 
class where g is substituted by 2g. The factor 2 is considered to be due to the chiral symmetry. 



C. Explicit results for a = 1 



The case a = 1 is specially interesting. On the one hand it is related to the spectral correlations at the Anderson 
transition and on the other hand the analytical results of the previous section are greatly simplified since the spectral 
determinant can be evaluated exactly. 

In the KM domain the summation in can be easily performed, 



^1 - 27r2 n2 ~ 12' 

n=l 



The DoS and the TLCF are consequently given by 



P2{ZI,Z2) 



1 



96.g2 
1 



2z ^ '^^ 

dz dz^ 



1 



4852 



d 



d 

^1^^ 1" ^2-^ 

ozi az2 



1 



9652 



,9^ 

az2 



d d 2 

2ziZ2 — h Z2 

ozi 0Z2 



dz^ 



K{zi,Z2) 



The results in AA's domain 1 <C zi_2, 1^5 are greatly simplified since both the spectral determinant 

z/25 



^^^-^ n 1 + {z/2TrgnY sinh 



(^/2ff)' 



(48) 



(49) 



(50) 



(51) 



and the momentum summations of the propagator can be exactly evaluated 

8^2g2 — _ iz/2Trgf ~ 



Re5:tf(g,e) = ^Re^ 

00 

Im y n(q, e) = — y , 

■■j^^ 4:11 g ^-^_^\n + iz 1271 g n — iz/2'Kg 



(l-V^z)) 



1 



1 

'2^ 



z z 
1 — — coth — 
2ff 2g 



(52) 



By utilizing the relations 



1 

2^ 



coth 



Zl + Z2 



coth 



Zl - Z2 



4zi 



43 4,9 
ViV2{V\VZ.^ - 1 



4ziZ2 
[zi ~ Z2) 



^{Vl~ViV2), 



(53) 
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we find the following simple form for the TLCF, 

sinfzi 



We easily see that the unitary limit zi - 
bracket. 



Z2, 



Z2 



cos(zi + Z2) 



V. 



(54) 



Zl — Z2 Zl+ Z2 

00, Ea. H24|) . can be recovered by keeping only the first term in the 



IV. CHIRAL SYMMETRY AND ANDERSON TRANSITION 



In this section, for the special value a = 1, we put forward an exact result for the DoS and the TLCF, valid for all 
z,zi^2, in the limit g ^ 1. The spectral fluctuations in this case present features as scale invariance, level repulsion, 
and asymptotically linear number variance similar to the ones at the Anderson transition. Thus the chRBM at a = 1 
is an ideal candidate to investigate the interplay between chiral symmetry and wavefunction localization leading to the 
Anderson transition. We recall that in the thermodynamic limit, according to the RG analysis of previous sections, 
the chiral symmetry will delocalize the eigenstates close to the origin and the DoS is divergent (though according to 
chiral random matrix prediction this divergence is cut off before the origin). 



A. Generalized chiral random matrix model for critical statistics 



We find exact expressions for the DoS and TLCF of the chRBM at a = 1 by mapping it onto an exactly solvable 
generalized chiral random matrix model (gchRMM) . The model in question is defined by the probability density 



P<(H) (X J VUexpTr 



1 

Ah 



where TC is given by Eq.(^, U has the chiral block structure 

' C/2 



(55) 



(56) 



with Ui,2 n X n unitary matrices, and DU denotes the Haar measure of the unitary matrices [/i,2- We shall see that 
the parameter h is indeed related to the conductance g hy h ^ 1/g. In Refs.|38l l39|. it was found that the spectral 
correlations of this model are equivalent to the spatial correlations of the Calogero-Sutherland (CS) model ^3 at 
finite temperature. By using this analogy, exact expressions for the spectral correlations were found for any g 1/h. 
It turns out that in the g ^ I limit such model is also equivalent to our chRBM at a = 1 The argument is as 
follows. By decomposing Ea. H55() into the blocks of H and U, the probability density can be written as 



P{n)oc J VUexpTr 



2^2 



1 + Sh-'n 

2/r 



-CC^ + 'in^hReUiCU2Cl 



(57) 



The integral over U is, in principle, performed by using the Harish-Chandra-Itzykson-Zuber formula. However one 
could attempt an alternate route. Random matrix theory predicts that the spectral density of U in the large n limit 
is constant and, due to level repulsion, the eigenvalues are consequently well separated. One can thus model the 
spectrum typical of a f7 as exp(i6'm) [m — 1, . . . 2n) where in a first approximation the phases 6m = 2'Km/2n are 
assumed to be equidistant. By doing that we are neglecting the small fluctuations around the equidistant position. 
Within these approximations the integral over U is reduced, up to multiplicative constant factors, to the contribution 
of one "typical" U . In a basis such that the Ui^2 in (|57|l becomes diagonal. 



P[H) (X exp i ^ IQ 



This expression describes a chRBM with variance 



1 + [Anhy sax 



2 „• 2 '^{i - j) 



2n 



1 + (4n/i)2 sin^ 1 + iir^h^r^ ' 



(58) 



(59) 
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but this is nothing but our original chRBM for a = 1 with rg ~ 1/h. Thus the spectral correlations of both models 
should coincide. Let us first review the exact results for the gchRMM reported in Ref . |35| . In the h <^ 1 limit, the 
DoS was found to be 



Ay_i/,. cosh^ [2t(i + ^)ri 'aJ' 



(60) 



where A — ixjln is the mean level spacing at ft, = 0, and p^°-'(z) is given by Ea. H31|) . The TLCF is 



((p(ei)pfe))) = - 



A2 



cosh^ [2t(l + f )] ^ A ' aJ\ ' 



where K{zi, z-i) is given by Ea. H45|) . 

The mean level spacing A(/i) at /i 7^ can be evaluated from Ea. H6U|l as 



(61) 



A(h) 



1 

A 



dt 



1/h 



cosh^ [2t{\ + ^)] A 



96 



(62) 



We note that this quantity corresponds to A in Ea. l29|) . Naive calculation of A using the nonlinear a model gives 
a divergent result and we thus find disagreement with Ea. (|62() . However, just like in Eas. (|30|l and H43() . if the above 
results are scaled in terms of A{h) as 



pi(z) = A{h){p{e = A{h)z/n)), 
P2{zi,Z2) = A2(/i)((p(ei ^A(ft)ziA)p(e2 = A(/i)z2A))), 



(63) 



we show below that agreement between both calculations is found. 

In order to prove this claim we perform a series expansion on the exact results (I6U|I and (|61|l to compare them 
with the findings of the previous section in both the AA and KM domain. For the KM domain we start with the 
expression (|6U|I . The integral is strongly peeked around t — 0. Therefore we can perform an expansion in powers of t 
up to terms involving corrections. 



dt- 



1 



1/h 



cosh^ 2t 



l + 2ht~ 2hr 



, sinh 2t 



dz 



1 



A{h) 



~96 



2z— + z' 



cosh 2t 

7T^ d. 

96 ^dz 



sinh 2t 
cosh 2t 



sinh^ 2t 
cosh^ 2t 



(0), 



dz 



dz"" 



(64) 



The TLCF is calculated in the same way. After a laborious calculation one is left with 



((P(ei)p(^2)» 



1 



+ - 



A2(/l) 
TT' 



48 



d 



zi 



Z2- 



96 



dzi dz2 

d d ry 

k 7.^ 

dzi dz2 



^ dz? ^ dz-i dzo ^ ^ dz? 



K{zi, Z2) 



(65) 



In the AA's domain 1 ^ z, zi 2, the DoS reduces to 



(P(e)> 



1 



A{h) 
1 



dt- 



1 



cosh^ 2t 



1 - 



cos [2z{l + ht)] 



2z 



cos2z„, , 1 



1 

8^2 



(66) 



where 



V(z) 



dt 



cos{2zht) TThz/2 
cosh^(2i) sinh(7r/iz/2) ' 



(67) 
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In a similar way, by using 

we obtain for the TLCF 

((P(ei)p(e2))> 



, . sm zi - Z2) cos(zi + Z2 
K(zx,Z2) ■ , (68) 

Z\ - Zl Zl+ Zl 



sin(zi - 2:2) ^ / Zl - Z2 \ _ cos(zi + Z2) ^ f zi + Z2^'''^ 



Z\ — Z2 V 2 / Zl + Z2 



(69) 



As expected these results are in complete agreement with the findings of the previous section by setting h — l/i^g. 
We note that the function 2?(z) in Ea. H23|l can be written in a spectral determinant form 

00 2 

^(^) - n ^^irkvi- ^''^ 

n—l ' 

We remark that the mapping proposed in this section provides with explicit exact results for any value of z. Con- 
cerning the relation with the spectral properties at the Anderson transition we mention that in R.ef.|35j it was shown 
analytically that the gchRMM (and consequently our chRBM at a = 1) reproduces all the features of critical statistics. 
We do not repeat here this discussion and refer to it for details and to the next section for numerical verification. 

Once we have proposed exact relations for the spectral correlations at the special value a = 1, it is worthwhile 
to ask whether an analogous result can be extended for the rest of as. Unfortunately a mapping as the one above 
reported for a = 1 cannot be extended to other values of a. We attempt a different strategy. The DoS l|()U|) and 
the TLCF (|61ll were explicitly evaluated in Ref.js^ in the large conductance g — l/nh '3> 1 limit. The resulting 
expressions, valid up to corrections, were derived keeping the combination z/g fixed. This amounts to neglect 
z^ I g^^"^ terms as compared to z^ j . We remark that this expansion coincides with the supersymmetry calculation 
in the AA domain but not in the KM domain where l/t;^ corrections were kept. The results of Ref . [ssjl (Eqs.(61) and 
(69) there) can be rewritten in the following way, 

Pi(^) = y(./o (^) + + \Uz).h{z){V{z) - 1), (71) 

2 



P2{x,y) = iRe^(n^+n2_^ ' ^+ 



2?i - 1 Pi - 1 



2(a; + y)2 2{x ~ yY 

2 



2(L^) ^+ [^i(^)-^o(y) + J^{x)JM\ + 2(f^y) ^- [M^)My) - M^)My)]^ , (72) 

where 'D± ='D{{x± y)/2) and 11 j_ = 11^(9, {x ± y)/2) with V and 11 given by Eqs.ljSlJ and lf?7|) at a = 1 respectively. 
These expressions are consistent with Eas. (|42|l and H54II . For a = 1, using the relation H52|l . the first line of the r.h.s of 
Ea. (|72|l vanishes. The reason of keeping the expression of the first line is that we want to separate purely perturbative 
contributions involving the propagator 11 from the nonperturbative ones. 

Since the dependence in a above is only through V and 11, we speculate that for general a the expressions for the 
DoS and the TLCF of our chRBM (in the g S> 1 limit) are given by the above expressions with T) and 11 replaced 
by Eas. H23(l and (|27|l respectively. For the DoS, this speculation is supported by the supersymmetric calculation 
(which is valid for arbitrary a) since Ea. (|71|l is consistent with Ea. l42|l . For the TLCF the situation is less clear. The 
supersymmetry result for the TLCF at a 7^ 1 has a complicated form 1461) due to the nontrivial mixing of perturbative 
and nonperturbative contributions. We remark that for the special value a = 1 it was greatly simplified to Ea. H54(l by 
using nontrivial relations of the propagator and the spectral determinant as demonstrated in Sec lIII CJl We could not 
manage to recast H4tifl as in Ea. H72l) . We thus conclude that Ea. (|72f) for arbitrary a is just the simplest alternative 
among several possible options. 



B. Numerical results 



We now compare the analytical predictions of previous sections with the results of direct numerical diagonalization 
of the chRBM in the spectral region close to the origin. Since we are mainly interested in the spectral correlations 
beyond the Thouless energy attention is focused on long range correlators as the number variance S^(i) = (i)^ — 

(L^) = J^^ dzpi{z) — J^^ dxdyR{x,y). We recall that the number variance S^(L) measures the stiffness of the 
spectrum. In the metallic regime, for eigenvalues separated less than the Thouless energy, fluctuations are small and 
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FIG. 1: DoS pi{z). The solid line is the analytical prediction of the chRBM (ITU at a given g. The symbols represent the 
numerical results for the indicated volume and bandwidth ro, and the thin line stands for the result at ergodic limit g ^ oo 
(chiral Gaussian unitary ensemble (chGUE)). As observed the agreement with the analytical results is impressive. 



S^(L) ~ log(i) for L 3> 1. Beyond the Thouless energy spectral fluctuations get stronger and ~ L'^^''^, where 

d is the dimensionality of the space. For disorder strong enough eigenvalues are uncorrelated (Poisson statistics) and 
E^(L) = L. At the Anderson transition, the number variance is asymptotically proportional to xL (x < 1) ^2l1|. The 
spectral fluctuations are studied by direct diagonalization of the chRBM for different sizes ranging from N ^ 2n = 240 
to iV = 2n = 1200. The eigenvalues thus obtained are unfolded with respect to the mean DoS. The number of different 
realizations of disorder is chosen such that for each N the total number of eigenvalue be at least 2 x 10^. In order to 
reduce unwanted finite size effects we have utilized a "periodic boundary" version of Eq. ^ as in Ref. |31| . 

We first study the spectral fluctuations at the critical value a = 1. In this case, since the spectral determinant 
can be expli citly evaluated, both the DoS (|71|l and the TLCF (|72|l have a simple form. It is straightforward to show 
analytically [33 that S2(L) - L/Anro for L » 1. In FigHwe show the numerical DoS versus the analytical prediction 
H71|l . As observed, the oscillations of the DoS are damped with respect to the ergodic regime. We stress that, since 
g — vq, the comparison between analytical and numerical is parameter free. 

One of the signatures of an Anderson transition is the independence of g on the system size. In our case, analytically 
it is also predicted that g — tq (ro ^ 1) is scale invariant. In Fig|2]the number variance for tq = 3 and a = 1 is plotted 
for different volumes N. As expected the spectrum is scale invariant as at the Anderson transition. In Fig|31we test 
our analytical findings by comparing them with the numerical number variance at the origin for different matrix sizes 
and bandwidth ro = 3. As observed, the theoretical predictions are fully confirmed. The deviations appearing after 
X ~ 20 are a well understood finite size effect. 

We now study the spectral correlations for non critical values of a. For 1/2 < a < 1, according to our conjecture 
Ea. (|72|l . the asymptotic value of the TLCF is controlled by the power-law tail of the perturbative part. Thus it 
is expected that for L ^ 1, ~ L^/^"~^. The dimensionless conductance g ~ rg"~^A^^~^" increases as the 

system size does. In Fig0|we show the numerical number variance for different as versus the conjecture Ea. H72|) . The 
numerical number variance follows the expected asymptotic power-law behavior in fair agreement with the conjectured 
result. This evidence is not conclusive since different conjectures with the same asymptotic limits also may be in good 
agreement with the numerical results. Clearly further work is needed to settle this issue. 
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FIG. 2: Number variance Y?{L) close to the origin of the spectrum. The symbols represent the numerical results for the 
indicated volume A'' and bandwidth ro = 2. As observed, the spectrum is scale invariant. Differences for large distances are 
due to finite size effects. 



V. APPLICATIONS: THE QCD VACUUM AS A DISORDERED MEDIUM 

In this section we discuss applications of our results in the context of QCD. Generally, for a > 3/2 corresponding to 
normal diffusion, our model should describe the leading finite g corrections to the spectral correlations of any disorder 
system with short range disorder and chiral symmetry. The case ofl/2<a<3/2 corresponds with chiral disordered 
systems with long range hopping. This is typical of systems driven by dipole interactions (l/r*^) 131 • We would also 
like to mention that for a = 1 our results provide a phenomenological description of the properties of systems with 
chiral symmetry at the Anderson transition. Before going into details a word of caution is in order. As mentioned 
in the introduction, the spectral correlations of chiral systems, unlike for the standard universality class, are highly 
dependent on microscopical details of the model as the exact form of the disordered potential. Thus the applications 
below reported must be considered as educated guesses among a broad class of systems in which our findings may be 
relevant. 

In the infrared limit the eigenvalue correlations of the QCD Dirac operator do not depend on the dynamical details of 
the QCD Lagrangian but only on the global symmetries of the QCD partition function Thus random matrices with 
the correct chiral symmetry of QCD (termed chiral random matrices) 6] accurately describe the spectral properties 
of the QCD Dirac operator up to some scale known as the Thouless energy. For larger energy differences, dynamic 
features become important and the standard random matrix model ceases in principle to be applicable. However, in 
a recent work ji^l , it has been reported that the spectral correlations of the QCD Dirac operator in a background of 
instantons are accurately described by a chiral RBM with a = 3/4 up to scales well beyond the Thouless energy. The 
reason for such exponent is indeed related to the dipole-like interaction dominating the interaction between quark zero 
modes and instantons. The matrix elements Tja (which physically describe the amplitude of probability for a quark 
to hop from an instanton to an antiinstanton) of the QCD Dirac operator in a basis of chiral zero modes decay as 
Tai ~ l/i?TA (dipole-like interaction) where Rix is the instanton-antiinstanton distance in a 4D space. It was shown 
in Ref. p4| that the spectral properties of systems with power-law hopping are similar in different dimensions provided 

that the decay exponent matches the space dimension. Thus a decay Tai ^ 1/^ia equivalent to Tai ~ 1/-^ia^ 

in ID. Furthermore, with this value of a, g scales at the bulk as a 4D conductor [g ^ y/n) which is the expected result 
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FIG. 3: Number variance T?{L) close to the origin of the spectrum for a = 1. The solid lines are the analytical predictions for 
the chRBM Eas. 11711 and 1721 . and the symbols represent the numerical results for the indicated bandwidth ro and A'^ = 400. 
As observed, the agreement with the analytical results is remarkable. The value of g chosen corresponds with the analytical 
prediction except for ro = 2 where the best fitting is for <; = 1.8 instead of g = 2. This deviation is expected since the analytical 
predictions are valid only in the <? ^ 1 limit. 



from chiral perturbation theory and lattice simulations j3| • The range of applicability of this chRBM was found 
not to be restricted to the above mentioned zero temperature case. As usual in field theory, temperature is introduced 
by compactifying one of the spatial dimensions. Thus in Euclidean QCD the effect of temperature is to reduce the 
effective dimensionality of the system to three. Now since the effective dimension of the space matches the power-law 
decay of the QCD Dirac operator ^/R^) one expects, according to Refs. p3.l45l |. multifractal wavefunctions typical 
of a metal-insulator transition. As mentioned previously, this situation corresponds with the case a = 1 in our chRBM. 
The above findings suggest that, in case that the restoration of the chiral symmetry expected at finite temperature 
were dominated by instantons, the physical mechanism leading to the quark-gluon plasma state of matter could be 
similar to an Anderson transition driven by dipole interactions. 



VI. CONCLUSION 



We have studied the spectral properties of a disordered system with chiral symmetry and long range hopping. By 
mapping the problem to a supersymmetric nonlinear a model we have obtained explicit expression for both the DoS 
and TLCF in different domains. It has been observed that, as in the nonchiral case, if the power-law decay matches 
the dimensionality of the space the spectral correlations are similar to the ones at the Anderson transition and are 
described by critical statistics. Based on the relation with a gchRMM we have put forward an exact formula for a 
specific value {a = 1) of the power-law decay and we have also speculated that a similar relation should hold for the 
rest of as. Finally we have argued that in the context of QCD, our model may be utilized to describe the spectral 
correlations of the QCD Dirac operator beyond the Thouless energy at zero and at finite temperature. 

A.M.G. was supported by the EU network "Mathematical aspects of quantum chaos". K.T was supported by 
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FIG. 4: Number variance S^(L) close to the origin of the spectrum in the delocalized regime a < 1. The symbols correspond 
with the numerical simulation. The lines correspond with the analytical number variance with R{x, y) given by the conjecture 
Eg. 1721 . In all cases A'' in set to TV = 400 and ro = 2. The value of g chosen is the best fit to the numerical results. We checked 
that it is close (±5%) to the analytical prediction, g ~ Tq^'^A^'^"^" [2^. 
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